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ABSTRACT 

Even though water is the main constituent in interstellar icy mantles, its 
chemical origin is not well understood. Three different formation routes have 
been proposed following hydrogenation of 0, O2, or O3, but experimental evi- 
dence is largely lacking. We present a solid state astrochemical laboratory study 
in which one of these routes is tested. For this purpose O2 ice is bombarded by 
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H- or D-atoms under ultra high vacuum conditions at astronomically relevant 
temperatures ranging from 12 to 28 K. The use of reflection absorption infrared 
spectroscopy (RAIRS) permits derivation of reaction rates and shows efficient 
formation of H 2 (D 2 0) with a rate that is surprisingly independent of temper- 
ed \ ature. This formation route converts 2 into H2O via H2O2 and is found to be 
orders of magnitude more efficient than previously assumed. It should therefore 
be considered as an important channel for interstellar water ice formation as 

O ' illustrated by astrochemical model calculations. 
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Subject headings: astrochemistry - - infrared: ISM - - ISM: atoms - - ISM: 
k> , molecules — methods: laboratory 

1. Introduction 

Solid water ice has been observed on the surfaces of many different astronomical objects. 
In the Solar System it is found on planets and minor bodies such as comets, trans-Neptunian 
objects and Centaurs. In dense, cold interstellar clouds, infrared observations s h ow that in 



terstellar dust grains are cove red with water-rich ices (jGillett fc Forrestl 1 19731 : iGibb et al. 



20041 ; iPontoppidan et al.ll2004L e.g.). The formation of these ice mantles is especially impor- 
tant in the process of star and planet formation, when a large fraction of heavy elements 
can be depleted onto grains. In the dense cloud phase water layers form on the bare grain 
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surfaces. Then during the gravitational (pre-) collapse, virtually all gas phase species freeze- 
out on top of these water layers resulting in a CO dominated layer that likely also contains 
traces of O2 but very little water. 

The observed H2O ice abundance cannot be explained by direct accretion from the gas 
phase only. The exact mechanism by which water ice is formed is not understood. The 
Herschel Space Observatory, to be launched in the near future, will provide important new 
information on gaseous water in interstellar space and will measure quantitatively the water 
abundance as a function of temperature, UV field and other parameters. Furthermore, the 
Photodetector Array Camera and Spectrometer (PACS) on Herschel will cover the 62 jum 
band of solid H2O. In this way Herschel will provide a unique opportunity to observe the 
bulk of the water bands that are unobservable from the ground and relate them to Spitzer 
and groundbased mid-IR observations of ices in protostellar envelopes and protoplanetary 
disks. Understanding the processes by which water forms and why it is not formed under 
other circumstances will be essential for the interpretation of these data. 



Tielens fc Hagenl (Il982l ) proposed a reaction scheme in which water ice is formed on the 



surfaces of grains via three different routes: hydrogenation of O, 2 , and ,^. Models predict 
that water can indeed be formed through such reactions in dense clouds (jTielens fc Hagen 



19821 ; Id'Hendecourt et al.lll985t lHasegawa fc Her 



DSt 



19931 ; ICuppen fc Herbst 



20071 . e.g.). Us- 



ing a Monte Carlo approach ICuppen fc Herbstl (120071 ) showed that the contributions of the 
different formation channels to water ice formation as well as its abundance strongly depend 
on the local en vironment. Howev e r, the initial reaction scheme with the corresponding rates 
as proposed by iTielens fc Hagenl (119821 ) is based on old, in some cases outdated, gas phase 
data of the equivalent reactions. Progress has been severely hampered by the lack of real- 
istic experimental simulations of these low-temperature, solid state reactions. Preliminary 
la boratory studies of w ater syn thesis testing the fi rst reaction channel have been reported 
by iHiraoka et al.l (119981 ) and by iDulieu et al.l (120071 ) . Both groups investigated the products 
of D- and O- reactions on an ice substrate (N2O and H2O, respectively) using temperature 
programmed desorption (TPD). In experiments exclusively using this technique, it is hard to 
rule out any H 2 formation during warm up. Furthermore, quantitative interpretation can 
be tricky because unstable species like H 2 2 are destroyed in the mass spectrometer upon 
ionization leading to an artificially enhanced H 2 0/H 2 02 ratio. 



The present work focuses on the H + O2 channel in which O2 is converted to H2O via 



H0O0: 



H + 2 — H0 2 

H + H0 2 - H 2 2 

H + H 2 2 ^H 2 + OH 



(1) 
(2) 
(3) 
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H + OH^H 2 



(4) 



According to ICuppen fc Herbstl (120071 ) this channel is, together with the 3 channel, re- 
sponsible for water formation in cold, dense clouds. The exposure of 2 ice to hydrogen 
and deuterium atoms is investigated by means of reflection absorption infrared spectroscopy 
(RAIRS) and TPD. These techniques allow one to determine formation yields and the corre- 
sponding reaction rates. The present work comprises a study of hydrogenation and deuter- 
ation reactions of 2 ice for di fferent temperatures b etween 12 and 28 K, i. e. , up to the 
desorption temperature of 2 ( lAcharyya et al.l 120071 ). The formation of H 2 and H 2 2 is 
observed at all temperatures. An optimum yield is found at 28 K. 



2. Experimental 

Experiments are performed using an ultra-high vacuum set-up (P < 5 x 10~ 10 mbar) 
which co mprises a m a in cha mber and an atomic line unit. The set-up is discussed in more 



detail in iFuchs et al.l (120071 ). The main chamber contains a gold coated copper substrate 
(2.5 x 2.5 cm 2 ) that is in thermal contact with the cold finger of a 12 K He cryostat. The 
temperature can be varied with 0.5 K precision between 12 and 300 K. A precision leak valve 
is used to deposit 2 (99.999% purity, Praxair) on the substrate. Ices are grown at 45° with 
a flow of 1 x 10~ 7 mbar s" 1 where 1.3 x 10~ 6 mbar s _1 corresponds to 1 Langmuir (L) s _1 . In 
order to compare results from different experiments, the thickness of the 2 ice is 75 L for all 
samples studied and the substrate temperature is kept at 15 K during the deposition. An 2 
ice of 75 L consists of roughly 30 monolayers. This thickness is chosen to exclude substrate 
induced effects. Because a diatomic homonuclear molecule like 2 is infrared in-active, gas 
phase 2 is monitored during the deposition by a quadrupole mass spectrometer (QMS). 
After deposition at 15 K the ice is slowly cooled down or heated (1 K min -1 ) until a selected 
temperature is reached. Systematic studies are performed for different temperatures between 
12 and 28 K. 



H(D)-atoms are pr oduced in a well characterized thermal-cracking device (ITschersich fc von Bonin 



19981 ; ITschersich! 120001 ) . A second precision leak valve is used to admit H 2 (D 2 ) molecules 
(99.8% purity, Praxair) into the gas cracking line. In each experiment the H + H 2 (D + D 2 ) 
flow through the capillary in the atomic line is 1 x 10~ 5 mbar s _1 and the temperature of 
the heated tungsten filament, which surrounds the gas cracking pipe, is about 2 300 K. The 
dissociation rate and the atomic flux depend on the pressure and temperature (ITschersich 
20001 ) and are kept constant during all the experiments. A nose-shaped quartz pipe is placed 
along the path of the atomic beam in order t o cool down H(D)-atoms t o room temperature 
before reaching the ice sample by collisions (jWalraven fc Silveral 119821 ). The H(D)-atomic 
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flux nearby the sample is estimated, within 50 %, as 5 x 10 13 cm -2 s" 1 . At temperatures of 
12 K and higher, no blocking of surface processes by the presence of H 2 is expected in the 
ice. 

The newly formed species after hydrogenation (deuteration) of O2 ice are monitored 
by RAIRS using a Fourier transform infrared spectrometer (FTIR) running at a spectral 
resolution of 4 cm -1 in the range between 4000 and 700 cm -1 (2.5 — 14 jum). Typically the 
ice is exposed to the H (or D) beam for 3 (or 2) hours and IR spectra are acquired every few 
minutes. 

Systematic control experiments have been performed in order (i) to unambiguously 
confirm that the products are formed by surface processes and not by gas phase reactions, 
(ii) to check that any water present in the system does not affect the final results and 
(Hi) to verify that water and H2O2 formation occurs in the solid phase after H(D)-atom 
bombardment and not by H2(D2) -molecules addition. For (i) co-deposition experiments 
are undertaken in which H and O2 are deposited simultaneously. Water is only formed if 
the surface temperature is below the desorption temperature of oxygen, confirming that the 
presence of the oxygen ice is required for this reaction sequence to occur. Point (ii) is verified 
by using inert initial substrates like N2 ice to estimate the background water contribution, 
as well as by using different isotopologues ( 18 C>2, 15 N 2 and D). Finally (Hi) is checked by 
using pure H 2 (D 2 )-beams, i. e. without any H(D) present. 



3. Results 

The formation of both H2O2 and H 2 ice is confirmed by the appearance of their infrared 
solid state spectral signatures. Figure [1] shows typical RAIRS results for hydrogenation and 
deuteration of O2 ice at 25 K. From top to bottom a time sequence of four spectra is plotted. 
These spectra are difference spectra with respect to the inital oxygen ice. However, since 
our initial o xygen ice only consists of 30 ML, no features due to the intrisically very weak 



O2 feature (jEhrenfreund et al.l Il992l ; iBennett fc Kaiser! 120051 ) are observed in the original 



spectrum. Both the H2O and H2O2 clearly grow in time as the H-fluence (H-flux x time) 
increases. Similar features appear for the deuteration experiment, although here clearly less 
D 2 is formed. After fitting the infrared spectra with a straight baseline the column density 
(molecules cm" 2 ) of the newly formed species in the ice is calculated from the integrated 



inten sity of the infrared bands using a modified Lambert-Beers equation (IBennett et al. 



20041 ) . In the range of our spectrometer, water ice has two candidate bands for determining 
its column density, at 3430 cm -1 and 1650 cm -1 (3 and 6 /im, respectively). Since the strong 
3430 cm" 1 feature overlaps with the 3250 cm -1 band of H2O2, the weak feature at 1650 cm -1 
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was chosen to quantify the water column density. Since literature values of transmission 
band strengths cannot be used in reflection measurements, an apparent absorption strength 
is obtained from a calibration experiment in w hich a pure water ic e layer desorbs at constant 



temperature until the sub-monolayer regime (lOberg et al.l 120071 ). The uncertainty in the 
band strengths remains within a factor of two. Quantification of H2O2 is done using the 
1350 cm -1 band. As it is experimentally very hard to deposit pure H2O2 ice, the apparent 
absorption strength has to be obtained indirectly by assuming the ratio of the integrated 
band strengths bet ween the two bands in transmittance to b e the same as in reflectance 



H 2 0/H 2 2 = 0.57 flGerakines et al.l Il995t iLoeffler et all 12009 ). The band modes of solid 



D 2 and D2O2 ices have systematic peak position shifts of ~400 cm -1 with respect to the 
H 2 and H2O2 bands. The column density for deuterated species is obtained in a similar way 
from a calibration experiment, while the absorption strength for D2O2 is estimated assuming 
that H 2 0/H 2 2 = D 2 0/D 2 2 . 

In Fig. [2] the column densities of water and H2O2 are plotted as a function of the H- 
fluence (atoms cm -2 ) for different substrate temperatures. Figure [3] gives the equivalent for 
the deuterated species. The results for H 2 2 and D 2 2 are found to be very reproducible 
(errors within the symbols). Due to their low column densities the relative errors for H 2 
and especially D2O are larger. The H2O2 and D2O2 results all show the same initial linear 
increase followed by a very sharp transition to a steady state column density, with the steady 
state value increasing with temperature. The water results show similar behavior although 
the transition is not as sharp and the temperature dependence of the steady-state value is 
not as clear. The observation that the results for all temperatures show the same initial 
slope means that the rate of the reaction to H2O2 is temperature independent. The final 
yield is however temperature dependent and this depends on the penetration depth of the 
hydrogen atoms into the O2 ices; at higher temperatures H-atoms can penetrate deeper. 
This is discussed in more detail later. 



During the preparation of this manuscript we received a preprint by iMiyauchi et al. 



(120081 ) who performed a similar experiment for one single temperature (10 K). Our results 
for 12 K turn out to be close to their results, apart from an absolute scaling due to different 
assumptions on the band strengths. 



4. Determining the reaction rates 

Reaction rates are obtained by fitting a set of differential equations to the time evolution 
curves of H 2 (D 2 0) and H2O2 (D2O2). Usually a diffusive mechanism is considered to 
construct these equations. For H2O2 that is formed from oxygen and converted into water 
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as given in Eqs. [1]|1] the rate equation would be 

^^± = (kt P + k% p )hn K (t)no 2 (t) 

~ (C P + C p 02 )^H(t)n H2 o 2 (t) (5) 

with k^ op and n x the hopping rate and column density of species X and k\ and k 3 the rate 
constants of reactions 1 and 3. 



A wad et al.l (120051 ) and iMiyauchi et al.l (120081 ) applied similar expressions to determine 



the rates of the formation of methanol and water, respectively. They solved these equations 
under the assumption that the atomic hydrogen abundance on the surface remains constant. 
In this way expressions for the surface abundances of the reaction products were obtained 
with only the effective reaction rates as fitting parameter 

™h 2 o 2 = ™o 2 (0)- — X —t (exp(-/3 3 £) - exp(-p\£)) , (6) 

Pi — P3 

with 

Pi = (C P + *&)fcin H . (7) 

Fitting this expression to the experimental data points in Figs. [2] and [3] gives a very 
poor agreement because the model has an exponential behavior whereas the experiments 
for several temperature values clearly do not. This would moreover result in a different 
rate for H + O2 for each temperature whereas the experimental curves show that the rate 
is independent of temperature. For this reason we decided to use a different model. We 
consider two regimes. In the first regime (t < t t ) the hydrogen atoms get trapped into 
the ice with an efficiency that is independent of temperature. Once an atom is trapped it 
can diffuse efficiently and find an oxygen molecule to react with. This results in a zeroth 
order rate. Once nearly all oxygen molecules within the penetration depth are converted to 
H2O2 (t > t t ), diffusion becomes rate limiting and a diffusive mechanism for the reaction to 
H2O2 and H2O is applied. Since now the ice is changed from an O2 ice to an H2O2 ice the 
penetration depth of the H-atoms into the ice will change as well. This model is described 
by the following set of equations 

n .PO2 ~ PH2O2 , PlPH 2 2 /-, / a ,» / \ 

"h 2 o 2 = Pit + - (1 - exp(-py)) (8) 

P0 2 P3 PO2 

m, = flg^fi- lzfgCzM ) (9) 

PO2 V P3 J 



for t < tt and 



P\ 

^H 2 2 = PH 2 2 Q , 1 a [ eX P(~^3(^ ~ t t )) ~ 6Xp(-^i(t - tt))] 

Pi _ P3 
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+ (P0 2 -PH 2 2 - nHaO a (*t))(l ~ exp(-/?i(t - £*))) 

+ n Ha o a (*t) (io) 

™H 2 = ™H 2 o0t) + (PH 2 2 - WHa0(*t)) X 

, #,exp(-ft(t - t t )) - frexp(-/? 3 (t - *t))* 



/?i-/3 3 



(11) 



for t > t t with px the penetration depth of H-atoms into ice X in units of column density, 
(3\ and (3[ the effective rates of reaction 1 and 03 the rate of reaction 3. j3\ and j3z represent 
effective diffusive rates that include both the diffusion rate and the reaction rate, whereas 
(3[ is mainly determined by the hydrogen flux times the efficiency of H trapping into the ice. 

Fitting this model to the data, the three rates are found to be independent of temper- 
ature. We therefore apply the same average rate to describe the results at all temperatures 
and only the two penetration depths are allowed to vary between the experiments. The re- 
sulting curves are indicated by the solid lines in Figs. [2] and |3j The obtained rates are given 
in Tabled! The penetration depths are the steady state values in Figs. [2] and [3J These clearly 
increase with temperature to very high values. This suggests that as the 2 ice reaches its 
desorption temperature (~30 K) the structure becomes more open and the O2 molecules 
more mobile, allowing the hydrogen atoms to penetrate deeply into the ice. The H2O2 struc- 
ture on the contrary is much more dense and rigid and consequently the H/D-atoms cannot 
penetrate more than a few monolayers even at the highest temperatures. The temperature 
effect is also much less prominent in cases where the ice is comfortably below its desorption 
temperature. 

A clear difference in penetration depth between the hydrogenation and deuteration 
experiments can be observed. For the oxygen penetration depth a difference of a factor of 
two is found whereas for hydrogen and deuterium peroxide the difference can be even as high 
as a factor of six. 



Like in the models by I A wad et al.l (120051 ) and iMiyauchi et al.l (120081 ). the rates 0\ and 



j3 3 are the products between the rates of the reactions and the hydrogen surface abundance. 
The latter is assumed to be constant during the experiment and is the overall result of 
accretion, desorption and reaction with both H and O2. A comparison of (3[ and f3\ indicates 
that the reaction with H only plays a minor role. 

The uncertainties in f3\ and $3 are mainly determined by the fit and are within 50%. The 
error in (3[ due to the fit is much less, but here the main uncertainty is determined by the layer 
thickness. The error in the calibration of the layer thickness is a factor of two. Assuming that 
every hydrogen atom that gets trapped reacts with O2 and considering that two hydrogen 
atoms are needed to convert O2 to H2O2, the trapping efficiency can be determined from f3[ 



and the flux (5 x 10 13 cm -2 s -1 ). This results in a trapping efficiency of ~10% for deuterium 
and ~20% for atomic hydrogen. This is very close to the stic king efficiency of hydrogen 
atom s to a water ice surface of ~30% under these circumstances (jAl-Halabi fc van Dishoeck 
20071 ). The high efficiency of reaction 1 is c onsist ent with r ecent theoretical studies of this 
reaction in the gas phase. Kie et al.l (120071 ) and Ku et al.l (120051 ) found that this reaction 
can proceed barrierless for certain incoming angles. The rates (3\ and (3% have very similar 
values. This suggests that also reaction 3 is very efficient and H 2 formation is only limited 
by the penetration depth of H 2 2 . 

The rate of reaction [ 3] show s no significant isotope effect. This is in contrast with the 
results by iMiyauchi et al.l (120081 ) who found a significant isotope effect for this reaction using 
the diffusive model to fit their 10 K results. If a large barrier is involved in reaction [31 barrier 
crossing would proceed via tunneling and an isotope effect is expected (Watanabe, private 
communications). The fact that we do not observe a (large) effect at 12-28 K either suggests 
that the barrier for this reaction is low or that other mechanisms for the formation of water 
should be taken into account in the model to fit the data. In a future paper we plan to 
address this question in more detail. For now the conclusion remains that the formation of 
water from 2 and H is very efficient. 



5. Astrophysical discussion and conclusion 

The hydrogenation and deuteration experiments presented in this letter show an efficient 
mechanism to convert 2 ice to H 2 2 and ultimately H 2 0. In the model that describes these 
experiments, the rate limiting steps for formation are the trapping of the hydrogen into the 
2 ice and the penetration depth of the hydrogen into the H 2 2 . Astrochemical models take 
reaction, diffusion and desorption barriers as input. Our data shows that the formation of 
at least H 2 2 proceeds via a reaction with a barrier that is much lower than the value of 
1200 K previously assumed by lTielens &: Hagenl (Il982f ). It should therefore be considered as 
a route for water formation on interstellar grain surfaces. 

Care should however be taken when extending these findings from a laboratory environ- 
ment directly to interstellar ices. The temperature independence of the reaction rate that is 
observed, for instance, is directly due to the unique property of the 2 ice that allows hydro- 
gen to penetrate its structure and thereby preventing desorption back into the gas phase. In 
interstellar clouds the mantle surfaces would not co nsist of pure Q 2 ice b ut of a mixture of 
different species with water as its main constituent (IWhittet et al.l Il998l ) . The structure of 
these "dirty" ices and the binding energies to it would govern the desorption and diffusion 
behavior of adsorbates like 2 and H. Experiments on the binding of atomic hydrogen on 
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water ice surfaces showed that the surface abundance is gener ally dependent on temperature 
( IPerets et al.l 120051 ; iHornekaer et al.l 120021 ; iDulieu et al.l 120051 ) . Also hydrogenation experi- 
ments of CO showed that the decreasing H coverage for increasing temperature results i n 
lower effective rates (IWatanabe et al.ll2003l 120061 ; iLinnartz et al.l 120071 ; ICuppen et al.l 120081 ) . 
Here we discuss the implications of the present work to two different interstellar environ- 
ments: hydrogenation of an apolar (water-poor) ice mantle after freeze-out and cold cloud 
conditions where ice is formed from direct deposition of H and O. 

Observations s how that interstellar ice mantles consist of polar (water-rich) and apolar 
(water-poor) layers iTielens et al.l (Il99ll ). Apolar ices are thought to form during freeze-out 
in the densest parts of the cloud. In dense gas, most of the atomic oxygen has converted to 
O2 which can become subsequently a constituent of this polar phase. The lack of observed 
H2O2 and H 2 means that most of the O2 in the ice does not react to produce H 2 02 or 
H2O. This laboratory work however shows that hydrogen atoms can penetrate deeply into 
O2 ices and then react. In apolar interstellar ices the penetration depth is not determined by 
oxygen ice but by the main constituent of the ice mantle, CO. We therefore have performed 
additional laboratory experiments in which a mixture of CO and O2 ice is exposed to H- 
atoms. These show indeed that only the top few layers are hydrogenated and that the main 
part of the ice stays intact in agreement with the observations. Details of these experiments 
are to be published in a future paper. 

In cold and translucent clouds H- and O-atoms deposit onto the grain simultaneously. 
The species can then react immediately and O2 is converted all the way to water. Here the 
penetration depth observed in the laboratory becomes unimportant. What can we learn from 
the present experiments then? The fast reaction kinetics justify treating H + O2 and H + 
HO2 in a similar way as H + H are treated in astrochemical models, with the difference that 
probably only hydrogen is mobile for reaction. The results further suggest that the continuing 
rea ctions leading to water proceed with high efficiency. The original grain surface network 
by (ITielens fc Hagenlll982l ) includes two more water formation routes: via OH and via ozone. 
Under dense cloud conditions the ozone route was proposed to be the most efficient as is 
shown in the left panel of Fig. HI This figure plots the contributions for the three different 
H 2 formation channels using a reaction network limited to the three water formation routes 
an d without any dissocia tion reactions. The model para meters are taken f rom m odel Ml 
by iRume fc Herbstl (120001 ) and the initial conditions from ITielens fc Hagenl (J1982J) are used 
(n(H2) = 5 x 10 4 , n(H) = 1.41 and n(O) = 2.36 cm -2 ). These calculations include a barrier 
for reactions 1 and 3. The present work shows that these barriers are negligible and the 
middle panel plots the same model calculation without these barriers. The figure clearly 
shows that the O2 channel has a major contribution to t he overall water formation rate. 
Recent laboratory experiments by ISivaraman et al.l (J20071 ) on the temperature-dependent 
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formation of ozone by irradiation of oxygen ices by high energy electrons suggest that mobile 
oxygen atoms prefer the O + O pathway over O + O2 even though the O/O2 ratio in the ice 
is very small. This implies the presence of a barrier for the formation of ozone. The right 
panel gives the model results when a small barrier of 500 K for this reaction is considered. 
The contribution of the ozone channel is now further reduced. 

Figure 0] further shows t he H?0? abu n dance , n(H 2 2 ), for all models as produced 



through the 2 + H route. iBoudin et al.l (119981 ) set an observational constraint of the 



H 2 2 ice abundance toward NGC7538:IRS9 of 5.2% with respect to the H 2 ice abundance. 
All models are consistent with this number. 
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Table 1. The reaction rates obtained by fitting the model. For the corresponding 

uncertainties see the text. 





01 
[molcc cm 2 s 1 ] 


Hi 

[s- 1 ] 


fa 

[ S -1] 


H + 2 
D + O2 


5.4 X 10 12 

2.5 X 10 12 


1.3 x 10-3 
3.3 x 10-3 


1.8 x 10"3 
2.7 x 10-3 



0.000 



8 910 




D,0 



4000 3500 3000 2500 2000 1500 

Wavenumber (cm ) 



1000 




3000 2500 2000 1500 
Wavenumber (cm ) 



1000 



Fig. 1. — RAIR spectral changes of the O2 ice as a function of H-atom (left) and D-atom 
(right) bombardment at 25 K. Spectra at a H(D)-atom fluence of (a) 4 x 10 15 , (b) 4 x 10 16 , 

are given. 
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Fig. 2. — The H2O2 (top) and H 2 (bottom) column densities as a function of time and 
H-atom fluence for different O2 ice temperatures. The symbols indicate experimental data, 
the solid lines represent the fitted model. 
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Fig. 3. — Similar to Fig. [2] for D2O2 (top) and D2O (bottom). 
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Fig. 4. — The contributions of the three different channels to H2O formation in three mod- 
els and the H2O 2 abun dance, n(H.202) from the O2 + H channel, (left) Old network by 
Tielens &: Hagenl (Il982l ). (middle) new network without barriers for reactions 1 and 3 and 
(right) new network with barrier for O + O2 reaction. See text for details. 



